Transition from weak to strong measurements by nonlinear quantum feedback control 



Jing Zhang, 1 ' 2 ' 3 -B Yu-xi Liu, 4 ' 2 Re-Bing Wu, 1 ' 2 Chun- Wen Li, 1 ' 2 and Tzyh-Jong Tarn 5 ' 2 

department of Automation, Tsinghua University, Beijing 100084, P- R- China 
2 Center for Quantum Information Science and Technology, 
Tsinghua National Laboratory for Information Science and Technology, Beijing 100084, P- R- China 
3 Department of Physics and National Center for Theoretical Sciences, 
National Cheng Kung University, Tainan 10101, Taiwan 
^Institute of Microelectronics, Tsinghua University, Beijing 100084, P- R- China 
5 Department of Electrical and Systems Engineering, 
Washington University, St. Louis, MO 63130, USA 
(Dated: May 10, 2010) 

We find that feedback control may induce "pseudo" nonlinear dynamics in a damped harmonic 
oscillator, whose centroid trajectory in the phase space behaves like a classical nonlinear system. 
Thus, similar to nonlinear amplifiers (e.g., rf-driven Josephson junctions), feedback control on the 
harmonic oscillator can induce nonlinear bifurcation, which can be used to amplify small signals 
and further to measure quantum states of qubits. Using the cavity-QED and the circuit QED 
systems as examples, we show how to apply our method to measure the states of two-level atoms 
and superconducting charge qubits. 
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I. INTRODUCTION 

Quantum feedback control [l| is one of the central parts 
of quantum control theory and applications [2j-[6| owing 
to its potential ability to improve the stability and ro- 
bustness of the system. Besides the extensive theoret- 
ical studies [7Hll| . recent rapid developments on sensi- 
tive measurements and manipulation techniques in atom- 
optical fL2T - [l4| and solid-state systems [l5[ have made it 
possible to implement quantum feedback control in lab- 
oratories. 

In atom-optical systems, atomic ensembles are put in 
high-Q optical cavities, so that the information of the 
states of the atomic ensembles can be extracted by the 
probe lights transmitted through the optical cavities. 
Photocurrents induced by the probe lights are processed 
by field programmable gate arrays to generate real-time 
control signals, which can be fed back to design the elec- 
tromagnetic fields imposed on the atomic ensembles [l2j . 
Possible applications of quantum feedback control in such 
systems include state stabilization [jfjj . entanglement 
production f^|, spin squeezing [l8j], and state discrim- 
ination fl3l Similar studies can be found in solid- 
state systems [T|j|, e.g., the superconducting circuit-QED 
systems, in which quantum feedback control has been 
proposed to cool and squeeze the motion of a nanome- 
chanical resonator |20l - |23j . In these studies, the position 
of the nanomechanical resonator can be measured by a 
single electron transistor [U or a rf-SQUID [22] |24j [2o| . 

Most existing theoretical studies on quantum feedback 
control are concentrated on linear quantum systems, i.e., 
the dynamical equations of the system are linear in the 



'Electronic address: jing-zhang@mail.tsinghua.edu.cn 



Heisenberg picture and the feedback controls are linear 
functions of the system states, which can often be re- 
duced to standard classical control problems, e.g., the 
Linear Quadratic Gaussian control problem [t| 101. How- 
ever, essential differences arise when we study nonlinear 
quantum systems: (i) linear quantum systems possess 
evenly spaced discrete energy spectra, while the distri- 
butions of energy levels in nonlinear systems are un- 
even or continuous; (ii) linear quantum system start- 
ing from a Gaussian state always stays in a Gaussian 
state during the evolution, while nonlinear quantum dy- 
namics generally distorts the wavepacket. Moreover, the 
presence of the relaxation and dephasing effects in non- 
linear quantum systems may give rise to inherent phe- 
nomena in nonlinear classical systems, e.g., chaos [26[ 
and bifurcation 12711 . These dissipation-induced nonlin- 
ear effects [111 |27|| have various applications in labora- 
tories. For example, nonlinear dynamical bifurcation of 
a rf-driven Josephson junction is used to amplify small 
signals [28|, and further applied to the readouts of the 
superconducting qubit states in experiments [HI, [29j] . 

Here we propose a method to mimic nonlinear dynam- 
ics using a harmonic oscillator with the feedback control. 
Such a proposal can be widely applied, especially in cir- 
cumstances where a nonlinear amplifier like the rf-driven 
Josephson junction is not achievable. It should be no- 
ticed that the manipulation of nonlinear effects via feed- 
back has been widely studied in classical systems [30|. 
However, to what extent can we apply this to quantum 
systems? In particular, an extensively discussed ques- 
tion is whether or not the nonlinear effects [3 11 ] can be 
produced in linear quantum systems by quantum feedback 
control. Here we address this question by examining a 
damped harmonic oscillator driven by a nonlinear feed- 
back control, and shows that using this "pseudo" nonlin- 
ear amplifier can read out quantum states of qubits. 
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Our proposal is motivated by the recent developments 
for the quantum state readouts in superconducting quan- 
tum circuits via a nonlinear amplifier [28| and the quan- 
tum information processing using the cavity QED effect 
(e.g., in atom-optical systems, cavity quantum-dot sys- 
tems, and the systems for the interaction between super- 
conducting qubits and the transmission line resonator). 
In contrast to the measurement of the quantum states 
using a nonlinear amplifier (e.g., to read out the states 
of a superconducting qubit using a rf-driven Josephson 
junction [Hj]), one merit for our study is that in the 
long-time limit we can analytically analyze the dynam- 
ics of the qubit-oscillator system even in the "nonlinear" 
regime, which makes it possible to see how to control the 
rate of the information extraction so as to balance be- 
tween the measurement sensitivity and the measurement- 
induced disturbance in the bifurcation readout regime. 
Our study can be applied to the measurement of the 
states of an atom inside a cavity, quantum dots inter- 
acting with single-mode cavity field, or superconduct- 
ing qubits in the circuit QED systems. Without loss of 
generality, below we choose the atom-cavity system and 
the circuit QED system as examples to demonstrate our 
method, because these two kinds of systems are more ex- 
perimentally controllable and developed very quickly in 
recent years. Thus our proposal might be more possible 
to be demonstrated in these systems. 

The paper is organized as follows: in Sec. [Ill we present 
our feedback control proposal by a general model in 
which a damped harmonic oscillator is measured by a 
homodyne detection and then driven by the output field 
of the feedback control loop. The model of a system, 
in which a harmonic oscillator is dispersively coupled 
to a qubit and driven by an outer feedback control cir- 
cuit, is discussed in Sec. IIII1 Also, the analysis of the 
bifurcation-induced qubit readout by a harmonic oscil- 
lator with "pseudo" nonlinear dynamics is given in this 
section. The applications to the atom-cavity systems and 
the circuit QED systems are shown in Sec. IIV1 and the 
conclusions and discussions are given in Sec. [Vj 

II. CONTROLLED DAMPED HARMONIC 
OSCILLATOR 

Consider a damped harmonic oscillator with an angu- 
lar frequency uj and a damping rate 7, which is driven by 
a control signal Ut . Let us assume that the decay of the 
harmonic oscillator can be detected by a homodyne mea- 
surement with efficiency 77. Under the weak measurement 
assumption [32[ , the dynamics of the harmonic oscillator 
with the state p and the measurement output dy for mea- 
suring the position operator x can be described by the 
stochastic equations [a-ll2|: 

dp = -UH,p]dt + ^[a]pdt + y/frfH[a]pdW, (1) 

dy = (x)dt + ^=dW, (2) 
V 2? 77 



where (x) — tr (xp) is the average of the position operator 
x; a and a* are the annihilation and creation operators 
of the driven harmonic oscillator whose Hamiltonian is 

H = hwa^a + hut x. (3) 

The term Hut x represents the interaction between the 
harmonic oscillator and a time-dependent classical con- 
trol field characterized by u t . The superoperators 2?[a]p 
and H[a]p are defined by: 

T>[a]p = a pa 1 — — (a) ap + pa 1 a) , 
H[a]p = ap + pa} — tr (ap + pa') p. 
The measurement is performed over the position operator 

x = — = (a + a 1 ) , 

whose conjugate momentum operator is 

p=-L(at-a). 

dW is a measurement-induced Wiener noise satisfying 

E{dW) = 0, {dWf = dt 

with £'(•) representing the ensemble average over the 
stochastic noise. 

As shown in Eq. ((2|), dW represents the measurement- 
induced noise in the measurement output of the homo- 
dyne detection. To reduce the influence of the noise in 
the feedback control design, we take the time average 

Yt =il!{ { * )dT+ ^ dw ) = il! dy w 

of the output signal, where t 3> I/7, as an estimation 
of the position x. Taking the long-time average of the 
stochastic signals is an effective filtering strategy to ex- 
tract stationary signals from the background noises, and 
the time-average signal, e.g., Y t in Eq. (@}, can be fur- 
ther used to design feedback control. Such a control de- 
sign has been used in the literature to prepare desired 
quantum states (see, e.g., Ref. [33| for the Dicke state 
preparation). 

With these considerations, we apply the following non- 
linear feedback control 

u t = f(Y t ) = -k 1 Y t + k 3 Y t 3 -k , (5) 

to the original system ^ described by the Hamiltonian 
in Eq. ([3]), where the positive numbers ko, k\, A; 3 are the 
control parameters. In the following discussions, the con- 
trol parameters k\ will be chosen such that k\ > 7. 

The nonlinear feedback control given by Eq. ([5]) in- 
duces a pitchfork static bifurcation by varying the bifur- 
cation parameter u near the bifurcation point: 

= ~ (h - • (6) 
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In fact, if the initial state of the harmonic oscillator is a 
Gaussian state po with 

(x) Po = 0, (p) Po = 0, 



(P 2 )po 



\ (pX + Xp) 



(P) 2 P0 = ^P0> 
( x )po \P)po = Ca; 



(7) 



term obtained by measuring a single system and used to 
update the feedback control. In contrast, our proposal is 
quite similar to a feedback control proposal based on the 
measurement over an ensemble of harmonic oscillators. 
In fact, the measurement output given by Eq. ((2|) can be 
reexpressed as: 



1(f) = (x) + 



1 



4(t), 



and if the bifurcation parameter uj is below uj*, then the 
state of the harmonic oscillator given by Eq. (JlJ converges 
to the following stationary coherent state 



Wo 



where 



2uj 

ry , 1 - 



UJ 2 - k\U) + 7 2 /4' 



(8) 



(9) 



However, if the bifurcation parameter uj exceeds oj*, the 
original stationary state given by Eq. ([5]) becomes unsta- 
ble, and two new branches of stationary coherent states 
appear 



, 00 \ 
h,2/ = 



1 



where 



2tu 
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+ k\0J - 7 2 /4 



k 3 uj 



(10) 



(11) 



The analyses of the above results can be found in Ap- 
pendix [X] It can be verified that jo^) are far away 
from the original stationary state |a§°) if the parameters 
ho, k 3 are small enough such that: 



ho, k 3 < 



uj 2 - k\0J + 7 2 /4 



(12) 



To give more insights about the above results, we see 
that Yt and x = E ((x)) coincide together in the long- 
time limit (see Eqs. (|A5|) and (|A8|) in Appendix |A")) , and 

thus Po°) and (xj° 2 > PT2) are just the stationary 
states of the dynamical equation: 
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x = --x 



up, 

p = — (oj — ki) X 



k 3 x 3 



(13) 



In fact, it can be verified that the above equation has one 
stable equilibrium (a;g°, p x ) T when u> < oj* and two sta- 
ble equilibria (xf 2 i PT2) when ui > uj* , which coincides 
with the results given in Eqs. © and (fTTj) . 

It should be pointed out that the feedback control pre- 
sented here is different from the Markovian feedback [8[ 
and Bayes feedback controls @ in which the white noise 
term dW in Eq. © is looked as an innovation information 



where £(£) satisfies 

=0, £(t)£(f) = <*(*-«')• 

Then, we take the average of the measurement output 
over the quantum ensemble to obtain x = E (I(t)) = 
E ((x)) which can be further used to design feedback con- 
trol ut = u (x) , and the same feedback control is imposed 
on each system in the quantum ensemble. Here, we have 
only a single system, and thus, different from the above 
ensemble feedback control proposal, we take the time- 
average Y t given in Eq. (Q| to replace the ensemble aver- 
age x. Since the system state given by Eq. ([1]) tends to a 
stationary state given by Eq. <j8j) or (fTO)) . Y t and x coin- 
cide together in the long-time limit from the ergodic the- 
ory [IJ (see Eqs. (|A5|) and (|A8|) in Appendix EJ. Thus, 
our feedback control proposal is similar to the ensemble 
nonlinear feedback control in the long-time limit. Com- 
pared with the existing feedback control proposals, e.g., 
the Bayes feedback, our method is more robust to the 
uncertainty of the initial state. In fact, the same station- 
ary states given by Eqs. (jU) and (fit))) can be obtained if 
the initial state of the system deviates slightly from the 
Gaussian state given by Eq. (J7J). 

To show the validity of our feedback control proposal, 
let us see some numerical examples. By setting the 
system parameters as: 7 = 250 MHz, kg — 50 MHz, 
k\ = 500 MHz, k 3 = 50 MHz, oj = 3.5 MHz, x = Pq = 0, 
V x „ = V Pa = 1, and C XaPo — 0, the simulation results in 
Fig. [IJi show that both (x) and Y t converge to a station- 
ary state ccg° = 0.0126 « given by Eq. ©. However, if 
we tune the oscillating frequency of the harmonic oscil- 
lator such that uj = 65 MHz, the stochastic trajectories 
of (x) and Y t are separated into two branches which tend 
to two different stationary states xf 2 — ±1-95 given by 
Eq. (fTTj) . The above simulation results coincide with the 
theoretical analysis. 

As analyzed in Appendix^ the first-order quadratures 
(x) and (p) of the controlled system evolve nonlinearly in 
the long-time limit, while the higher-order quadratures 
approach to those of the linear quantum systems which 
preserve the Gaussian properties of the states. Therefore, 
the nonlinear dynamics induced by the proposed quan- 
tum feedback control is "semiclassical" in some sense (see 
the simplified diagram of the feedback control circuit in 
Fig. [5]), in contrast to the dynamics of the system gov- 
erned by a fully quantum nonlinear Hamiltonian (about 
x and p) such as: 



nl 



1 2 1 1 2 2 ; 4 
— V H OJ X — kx 

2 F 2 



(14) 
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FIG. 1: (Color online) Stochastic pitchfork bifurcation given 
by Eqs. ©, ©, and © with (a) uj = 3.5 MHz and (b) cj = 65 
MHz. The parameter r = 2 ns is the normalization time scale. 
The red curves with plus signs (one curve in Fig. []J, and four 
curves in Fig. [lb) represent the stochastic trajectories {x)t, 
and the blue solid curves (one curve in Fig. [T^ and four curves 
in Fig. [Tp) denote the time-average trajectories Yt given by 
Eq. @. 
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Local Oscillator 



FIG. 2: (Color online) Schematic diagram of the feedback 
control circuit: the harmonic oscillator is measured by a ho- 
modyne detection, and the output signal of the homodyne 
measurement is fed into a classical control circuit which is 
composed of a integral filter and a controller. The output 
signal of the control circuit is further fed back to drive the 
harmonic oscillator, and this feedback control signal may in- 
duce nonlinear dynamics of the harmonic oscillator. 



in which the Gaussian wavepackets are distorted. Here 
the subscript "nl" means that the Hamiltonian contains 
higher-order nonlinear terms. In our feedback control 
proposal, we introduce nonlinear terms like (x) 3 to the 
dynamical equation of the system in comparison to the 
terms (x 3 ) introduced by the nonlinear Hamiltonian H n \ 



in Eq. (| 14|) . Thus, under the nonlinear Hamiltonian H n \, 
the equations of the first order quadratures (x) , (p) and 
the second order quadratures V x , V p and 



C xp are not 



closed, since higher-order quadratures such as (x ) are 
involved. However, in our proposal, the equations of (a;), 



{p} and V x , V p , C xp are closed and decoupled just like 



linear harmonic oscillator if the initial state of the har- 
monic oscillator is a Gaussian state (see the analysis in 
the Appendix [X]) , whose stationary solution can be an- 
alytically solved. For this reason, we call this feedback- 
control-induced nonlinear dynamics as "pseudo" nonlin- 
ear dynamics. As shown below, interesting phenomena 
can be observed when such a "pseudo" nonlinear system 
is coupled to another quantum system. That is, similar 
to the nonlinear amplifier using the nonlinear d yna mical 
bifurcation (e.g., rf-driven Josephson junctions [28j), the 
harmonic oscillator with feedback control can be used to 
amplify small signals, and furthermore read out the qubit 
states. 



III. MEASUREMENTS ON QUBIT STATES 
USING CONTROLLED HARMONIC 
OSCILLATOR 

A. Qubit-oscillator coupling 

To demonstrate how we can use the "pseudo" nonlin- 
ear harmonic oscillator produced by feedback control to 
measure qubit states, let us consider a coupled qubit- 
oscillator system driven by a classical field with angular 
frequency uid- The Hamiltonian of this system can be 



expressed as: 



hljj n 



H = ~~^~ a z + hiij a i a + hg(a) (j- + aer + ) 



y/2 



ae 



iui d t 



(15) 



where ui q and lu are the angular frequencies of the qubit 
and the harmonic oscillator; a z and a± are the z-axis 
Pauli operator and ladder operators of the qubit; a and 
cl* are the annihilation and creation operators of the 
harmonic oscillator; g is the qubit-oscillator coupling 
strength; and ut is the coupling constant between the 
driving field and the harmonic oscillator. 

In the dispersive regime, i.e., |A go | = \ui q — cj | 3> \g\, 
we can apply the following unitary transformation 



U = exp 



9 

A q0 



to the Hamiltonian H in Eq. (jT5j) , then we have an effec- 
tive Hamiltonian 



H = UHU^ 



h ( w 9 + X) 



a z + hui a'a + hxa aa 



-hu t 



V2 



(16) 



where \ — 9 2 /A qo is the effective coupling strength be- 
tween the qubit and the harmonic oscillator. 

Furthermore, the Hamiltonian H can be rewritten in 
the rotating reference frame under the unitary transfor- 
mation 

[/rot = exp [i (u d t) a* a] , 
thus we can obtain the following effective Hamiltonian 



H, 



c fr 



u Iot Hul t + iU Iot a 



rot u rot 



a z + hAoda^ a + hx<^ aa z + hu t x, (17) 
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where A d = uj q — cud is the angular frequency detuning 
between the harmonic oscillator and the driving field; and 
x = (a + a^) I \/2 is the position operator of the harmonic 
oscillator. Here the frequency-shift of the qubit caused 
by the harmonic oscillator has been neglected under the 
assumption that u) q 3> \ which is usually valid in atom- 
optical systems and superconducting circuits. 

If we expand the Hamiltonian in Eq. (|TB)) to (g/A qo ) 3 
terms, we can obtain the following effective nonlinear 
Hamiltonain: 




The last second term in Eq. (I18[) may induce nonlinear 
dynamics to the harmonic oscillator. However, (g/A qo ) 2 
times smaller compared with the third term in Eq. (|18[) , it 
is too small to be observed in the large-detuning regime 
g <C A qo . Thus, the high-order terms have been ne- 
glected in below discussions. 

The decay of the harmonic oscillator is detected by 
a homodyne measurement. Thus, the evolution of the 
qubit-oscillator system is conditioned on the measure- 
ment output of the homodyne detection, which can be 
described by the stochastic master equation ([T]) by re- 
placing the system Hamiltonian H with H e g given in 
Eq. (fT"7]) . The output signal of the homodyne detection 
is integrated to obtain a new output signal Y t — J Q dy/t, 
which is further fed into a nonlinear controller to pro- 
duce the control signal Ut given in Eq. ([5]) . As presented 
in Sec. UH such a simple third-order nonlinear feedback 
control induces a static bifurcation that can be used to 
enhance the measurement strength for the qubit readout. 

Different from the open-loop control that is predeter- 
mined by the designer without any information extrac- 
tion, the proposed feedback control u t is automatically 
adjusted according to the state of the qubit in real time so 
that different controls generate different state trajectories 
and thus different output signals. This makes it possible 
to identify the state of the qubit by amplifying the differ- 
ence between the output signals by the designed control. 
This feature of feedback control has been reported in the 
literature to enhance the measurement intensity by lin- 
ear amplification (see, e.g., Ref. (35[), which can be done 
more efficiently via nonlinear amplification induced by 
the proposed quantum feedback control. 

Let us assume that the qubit-oscillator system is ini- 
tially in a separable state 

p(0) = p fl (0)®|Vk>(0)>M(0)|, 

where 

P,(0)= E Pij(°)\i)(j\ 

is the initial state of the qubit with ground state \g) 
and excited state \e); and |-0o(O)) is a Gaussian state of 



the harmonic oscillator with the first and second-order 
quadratures given in Eq. (JTJ) . Then, the stationary state 
of the qubit-oscillator system in the long-time limit can 
be expressed as [36j : 

p°°= E p«wi ® i«?°><«n ( i9 ) 

where both and |a^°) are coherent states. Let 

oo ^ / ' oo _|_ ■ oo \ 

then [xf,pf) and (xf,pf) are respectively the sta- 
tionary states of the equations: 

T 

Xg,e = --X g , e + (A od T x)Pg,e, (20) 
Pg.e = ~^Pg. e - (A od TX~ k^Xg^e - k 3 X g e + k . 

The coefficients pf- in Eq. (fT9j) are given by: 

p™= Pgg (Q), P Z=Pee(0), p~=p~ = 0. 

As analyzed in the last paragraph of Sec. UH the har- 
monic oscillator can be looked as a "semiclassical" pseudo 
nonlinear system driven by an outer classical feedback 
control circuit (see the simplified version of the feedback 
control circuit in Fig. [3]). The interaction between such a 
pseudo nonlinear system and the qubit brings two aspects 
of effects. On the one hand, for the qubit, this interac- 
tion brings additional decoherence. In fact, it can be 
found that there exists a measurement-induced dephas- 
ing factor for the reduced states of the qubit of which 
the damping rate can be approximately estimated in the 
long-time limit as (see the analysis in Appendix [Bj : 

T = X (xTpT-pTxT). (21) 

On the other hand, this interaction leads to additional 
frequency shift for the harmonic oscillator depending on 
the state of the qubit. Therefore it is possible to disper- 
sively read out the states of the qubit under appropriate 
conditions. 



B. Bifurcation-induced quantum measurement 

The equations of the harmonic oscillator corresponding 
to the ground and excited states of the qubit given in 
Eq. (|20|) have the same form as in Eq. (IT3|) with two 
different angular frequencies u = A 0( j — \ an d w = A Q( j + 
X- As analyzed in Sec. |Hl a static pitchfork bifurcation 
occurs at a critical angular frequency 

= \ (fci - v^i-r 2 ) • 

When u < w*, the system given by Eq. (fi"3| possesses 
only one stable equilibrium, which becomes unstable 
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FIG. 3: (Color online) Simplified version for the feedback 
control circuit of the coupled qubit-oscillator system. The 
qubit is coupled to the harmonic oscillator, which is driven 
by a feedback control circuit composed of a filter and a con- 
troller. The nonlinear feedback control coming from the con- 
trol circuit leads to "semiclassical" nonlinear dynamics of the 
harmonic oscillator. The static bifurcation induced by such 
nonlinear dynamics can be used to amplify small signals and 
thus measure the states of the qubit. 



when u) > lu* and, in the meanwhile, two other stable 
equilibria appear. 

Therefore, by tuning the angular frequency uid of the 
driving field near the bifurcation point, the static bifur- 
cation introduced by the proposed feedback control can 
induce a transition from weak to strong measurements 
for the qubit readout. Actually, if 



uj d = uj - u* + 2x, 



(22) 



both effective angular frequencies uu g = lu* — 3% and uj e — 
lu* — x corresponding to the two eigenstates of the qubit 
are lower than lu*. Then, in the long-time limit, e.g., 
when t S> I/7, the measurement outputs corresponding 
to the two eigenstates of the qubit are 



Y g ,e(t) 



W L ^ kiuj g . e + 7 2 /4 



(23) 



respectively (see Eq. ([9]) and the analysis in Appendix 
rA"|) . If ko is small enough such that 



k < 



X 2 - 2 X oj* + k lX 



(24) 



both Y g (t) and Y e (t) are so close to zero that they are 
almost indistinguishable. When i > 1/7, we have 



Pg,e -> Pg, 



9,e 



lu 2 £ - fciw ff , e + 7 2 / 4 ' 

1 

2uu 9,e ' 



then we can calculate the measurement-induce dephasing 



rate r weak from Eq. (|2lj): 



weak 



(25) 



The tiny separation between Y g (t) and Y e (i) can be 
amplified by switching the angular frequency of the driv- 
ing field from LUd — u a — lu* + 2x to Cud — lu — <*)*■ Un- 
der this condition, the two effective angular frequencies 
ujg = oj* — x an d w e = lu* + x are lower and higher than 
the bifurcation angular frequency lu* respectively. When 
t I/7, the corresponding stationary measurement out- 
puts are 



Ygit) -> /,„- h . 



Ye(t) -► X% 

By setting 

k 3 < 

we have 



'0)2-^+72/4' 

-Q 2 + kxCje - 7 2 /4 
fc 3 w e 



X 2 +2 X lu* - hx 



lu* + x 



(26) 



(27) 



\x?-x?\*>\x?-x?\. 

The above result means that the difference between the 
two output signals corresponding to the two eigenstates 
of the qubit is amplified by the proposed nonlinear feed- 
back control. When t > 1/7, it can be calculated that 



Pb 

Pe 



x°° 

Xg 



pT 



= k - 



Cu 2 -k^g+^/A 1 



Cu 2 e + fcic^ e - 7 2 /4 
k 3 Cu e 




Then, we can obtain the measurement-induced dephasing 
rate r str0 ng from Eq. (l2lT) as: 



strong 



= xKpT-pT^T) 



fco7X 2 _ V-^e + kiu e - 7 2 /4 



(28) 



It can be noticed that r strong is far greater than Tweak 
when ko and £3 are sufficiently small to satisfy Eqs. (f2"4")l 
and IJ57J}. 
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IV. APPLICATIONS 



to be 



A. Atom-optical systems 

Consider a coupled atom-cavity system, in which a 
two-level C s atom is dispersively coupled to a single-mode 
field in an optical cavity [39[ (see Fig. 3]). The probe light 
transmitted through the optical cavity is detected by a 
homodyne measurement. Then, the output photoncur- 
rent is fed into a control circuit composed of an integral 
filter and a controller to generate a control signal, which 
is further fed back to control the probe light by a modu- 
lator. 

Controller Filter 









\_ 






Modulator 


— D/A — 


«,=f(X) 




— AID — | 



Probe 
laser 



Detector 




Mirror 



Local Oscillator 

FIG. 4: (Color online) Schematic diagram of the atom-optical 
system under feedback control. A two-level C s atom is cou- 
pled to a single-mode quantum field in the optical cavity, 
whose output field is detected by a homodyne measurement 
("LO" denotes a local oscillator). The output signal of the 
homodyne measurement is fed into an electric control circuit 
to generate the desired control signal, which is further fed 
back to modulate the probe laser fed into the optical cavity. 
The output signal of the homodyne detection may also be 
converted into a digital signal via an analog/digital (A/D) 
signal convertor, and thus the control signal can be gener- 
ated by a Digital Signal Processor (DSP), which is further 
converted into an electric signal via a digital/analog (D/A) 
signal convertor and fed back. 



k /2iT = 1 MHz, 
fci/27r = 6 MHz, 
fc 3 /27r = 1 MHz. 

We first tune the detuning A a d = uj q — uid between the 
frequency uj q of the optical cavity and the frequency ujd 
of the driving field such that 

A od /27r = -3.57 MHz. 

In this case, the long-time limits of the measurement out- 
puts Y g and Y e corresponding to the ground and excited 
states \g) and |e) of the two-level atom are both close 
to zero, and thus almost indistinguishable (see Fig. [5^). 
This corresponds to the weak measurement case. 
If we tune the detuning frequency A oc i such that 

A od /27r = 0.083 MHz, 

the two branches of outputs corresponding to \g) and |e) 
are separated in the long-time limit, which corresponds 




O -2 



(see Fij; 


;.Eb). 




Y e (tMe> 




r ' " 






(b) 





t/T 



2.5 5 7.5 10 



FIG. 5: (Color online) Measurement outputs of the atom- 
cavity system under feedback control, (a) A d/2n — —3.57 
MHz (weak measurement case): the two output signals Y g (t) 
(red curve with plus signs) and Y e (t) (blue solid curve) cor- 
responding to the two eigenstates \g) and |e) of the two-level 
atom are almost indistinguishable; (b) A d/27r = 0.083 MHz 
(strong measurement case): the two output signals Y g {t) (red 
curve with plus signs) and Y e (t) (blue solid curve) are sepa- 
rated in the long time limit. The parameter r — 0.16 fis is 
the normalization time scale. 



We choose the followin g expe rimentally accessible pa- 
rameters (see, e.g., Refs. |39l44l1 |): 

A go /27r = 35 MHz, 
g/2n = 8 MHz, 
7/2tt = 1.4 MHz, 

where A qo = tu q — lu is the detuning between the fre- 
quency uj q of the two-level atom and the frequency uj 
of the single-mode field in the optical cavity; g is the 
coupling constant between the two-level atom and the 
optical cavity; and 7 is the damping rate of the optical 
cavity. The control parameters fco, fci, and are chosen 



B. Superconducting circuits 

The second example considers the superconducting 
circuit shown in Fig. [6l in which a transmission line 
resonator is capacitively coupled to a Cooper pair box 
(charge qubit). Let us now discuss experimental feasibil- 
ity via the numerical simulation using the experimentally 
accessible parameters. According to the current experi- 
ments (see, e.g., Ref. [42j]), the qubit frequency u> q , the 
frequency lo and the decay rate 7 of the transmission line 
resonator, as well as the coupling constant g between the 
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Output 



-k 




Charge qubit 



LO 



FIG. 6: (Color online) Schematic diagram of the supercon- 
ducting circuit under feedback control. A charge qubit is 
capacitively coupled to a transmission line resonator (TLR), 
whose output field is detected by a homodyne measurement 
("LO" denotes a local oscillator). The output signal of the 
homodyne measurement is fed into an electric control circuit 
to generate the desired control signal, which is further used 
to drive the electric field in the transmission line resonator. 



qubit and the transmission line resonator can be chosen 
as: 



uiq/2ir 

7/27T = 



= 5.1 GHz, 
100 MHz, 



w /2tt = 5 GHz, 
g/2iT = 20 MHz. 



(29) 



With the conditions given in Eqs. (j2"4"]l and (|2"T|) . we 
further assume that the parameters ko, hi, and fc 3 are 



fc /27r = 20 MHz, 

fci/2vr = 200 MHz, 

fc 3 /27r = 2 MHz or 10 MHz. 



(30) 



The frequency uid/^n of the driving field is initially cho- 
sen to be 4.995 GHz. At the time t* /r = 100, i.e., t* = 50 
ns, the frequency u>d/2ir of the driving field is switched 
from 4.995 GHz to 4.987 GHz, where r = 0.5 ns is a nor- 
malization time scale. Simulation results in Fig. [7] show 
that, at time t*, the difference between the two output 
signals Y g (t) and Y e (t), i.e., the measurement sensitiv- 
ity, suddenly jumps. The measurement-induced dephas- 
ing rate also suddenly jumps from r wcak to r strong . In 



fact, it can be calculated from Eqs. (|25|) and (|28| that 
r wC ak/27r w 0.36 MHz and r stro „g/27r « 10.22 MHz (or 
4.57 MHz) when fc 3 /27r = 2 MHz (or 10 MHz). This in- 
dicates that the static bifurcation introduced by our pro- 
posal induces a transition from weak to strong measure- 
ments at time t*. Moreover, as shown in Eq. (|2"6")l and the 
simulation results in Fig. the decrease of the nonlinear 
coefficient k 3 makes the measurement more sensitive, but 
accelerates the dephasing of the qubit. Such a tradeoff be- 
tween measurement sensitivity and measurement-induced 



dephasing effects is a natural consequence of the conflic- 
tion between information extraction and measurement- 
induced disturbance, which is inherent for quantum mea- 
surement. 



gi 1 
to 
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FIG. 7: (Color online) Bifurcation-induced transition from 
a weak measurement to a strong measurement with (a) the 
two output signals Y g (t) and Y e (t) corresponding to the 
two eigenstates \g) and |e) of the charge qubit, and (b) 
the measurement-induced dephasing factor. The parame- 
ter t = 0.5 ns is the normalization time scale. The solid 
curves represent the trajectories with a nonlinear coefficient 
kz/liv = 2 MHz, while the solid curves with plus signs repre- 
sent the trajectories with kzj2-K — 10 MHz. 



V. CONCLUSION 

In summary, we present a method to induce "semiclas- 
sical" nonlinear dynamics in a damped harmonic oscilla- 
tor via quantum feedback control. The nonlinear feed- 
back control induces a static bifurcation, which can be 
used to amplify the small frequency-shift of the harmonic 
oscillator for the qubit readout. Theoretical analysis and 
numerical simulations show an evident transition from 
weak to strong measurements near the bifurcation point. 
Our proposal works as well as the bifurcation readout 
proposal by a nonlinear amplifier, e.g., a rf-driven Joseph- 
son junction. Additionally, we can tune the information 
extraction rate to balance between information extrac- 
tion and measurement-induced disturbance in the "non- 
linear" readout regime. Using the atom-optical systems 
and circuit- QED systems as examples, we show how to 
apply our proposal to experimental systems. 
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We emphasize that the proposed cavity-assistant non- 
linear amplification strategy can also be applied to the 
qubit readouts in other experimental implementations, 
such as quantum dot-cavity systems. This generalized 
result is important, in that the proposed measurement 
method might be more efficient, because the nonlinear 
amplification device, like the rf-driven Josephson junc- 
tions in superconducting circuits, may not be achievable. 

Our study is also hopeful to be extended to quan- 
tum demolition measurements, e.g., the detection of mo- 
mentum or position of a nanomechanical resonator, for 
which the measurement-induced back-action effects on 
the quantity being measured cannot be neglected. 
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are the variances of the position x and momentum p; and 

. xp + px . 



*><P) 



is the symmetric covariance of x and p. 

Taking the average of the above equations, we have 

x = -~x + up, 

- 7 - - s 

p = -ujx - -p + kiY t - k 3 Y t + ko, 

Y t = ~(Y t -x), (A3) 

where x — E ((x)), p= E ((p)), and Y t = E (Y t ). Here, 
we neglect the high-order correlation of Y t induced by the 
classical noise dW which is reasonable when the evolution 
time is sufficiently long. 
When w < w*, we have 



cu 2 - kiu) + 7 2 /4 
7/2 

oj 2 - kiui + 7 2 /4 



>0, 
> 0. 



(A4) 



Appendix A: Dynamics of the Controlled Harmonic 
Oscillator 

From the stochastic master equation (flj and the mea- 
surement output equation @ given in Sec. [Til we can 
obtain the following non-Markovian stochastic integro- 
differential equations for (x) and (p): 



Then, since we start from the initial state (0, 0, 0) , it 
can be verified from Eq. (|A3[) that 

0<x<x^, 0<p<p™, Q<Y t <x^. 



Thus, from the last equation in Eq. (| A3|) . we have Y t — > 0, 
i.e., 3Y°°, s.t., Y t -> Y°° when t -> oo. Substituting Y°° 
into the first two equations of Eq. (| A3|) , we have 



d(x) 
dip) 



^(X)dt+L0(p)dt+ \V X - i 

7, 



-uj{x)dt - j;(p)dt + k Q dt + y / 2i]"fC xp dW 



~ki - J |(.r}<// 



1 



-*3 



)dt 



1 



-.dW 



.dW 



n 3 



(Al) 



which are equivalent to stochastic Markovian equations 
by introducing a new variable Y t : 



d(x) 
dip) 

dY t 
where 



-l{x)dt + u(p)dt + (v x - M dW, 

-ui{x)dt - ~(p)dt + (hiY t - k 3 Y t 3 + fc ) dt 
-y/*f/C xp dW, 
l(Y t -(x))dt-j-±=dW, 



(A2) 



Since the matrix 



-u ~l J \ p 




< 0, 



we know that there exists (x 00 ,^ 00 ) such that 



(x{t) 
\P{t) 



V 



when t — ^ oo. 

Next, we will show that Y°° = x 00 . Note that 



Y t = E(Y t ) = - / x{r)dT, 



we know that Ve > 0, 3T a , s.t. Vt > T a , \x - x°°\ < e/2. 
Let 



V x = (x 2 ) - (x) 2 , V p = (p 2 ) - (pf 



MT a ) = max \x - x c 

[o,r„] 
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and 

TV =max{T (T ,2e- 1 T <T A(T a )} 
Then, Vf > Ty, we have 

|?i -ar°°| < - f \x{s)-x°°\ds 



a stationary coherent state ) given by Eq. ([8]) when 
oj < u> . 

With the same discussions, it can be verified that there 
exist two stable equilibria for Eq. (| A3|) when u> > oj* 
which can be expressed as: 



1 



pi ^ r 2+ ^: 7V4 i m 7| , ( A8) 



1 



|i( a ) 



Ids 



< 



which means that Y c 



In order to solve (a; 00 ,^ 00 ) , we let x — 0, p = to 
obtain the algebraic equations 



= 
= 



-x 



7, 



- top" 



+ k lX °°-k 3 {x°°) 3 + k . 



and Eq. (IA7|) can also be obtained in this case. It means 
that the state of the controlled harmonic oscillator tends 
to \af) or \af) when oj > oj*. 



Appendix B: Approximate estimation of the 
measurement-induced dephasing rate 

With the same analysis as in Appendix A of Ref. 43] , 
the state of the qubit-oscillator system can be expressed 
as: 

p(t)= E M*)i*>oi®hM*)x^(*)i> 



Although three equilibria can be solved in this case, only 
one equilibrium (a;g° is stable when oj < oj* , where 

are given by Eq. (|A4j) . Thus, if oj < oj*, we have 



E((x)) 
E((p)) 
E(Y t ) 



P0; 

+ Y n c 



b ' 



(A5) 



when t — v oo . 

Furthermore, it can be calculated that the equations of 
higher-order quadratures are independent of the control 
Ut- Thus, the evolutions of higher-order quadratures are 
like those of linear harmonic oscillator. Since we start 
from a Gaussian state, the state given by Eq. (p} remains 
to be a Gaussian state, and thus we only need to calculate 
the variances V x , V p and the symmetric covariance C xp . 
From Eq. (fTJ, we can obtain that 
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-jV x 
-7V 



7 ( 1 

2ojC xp + - - 2?77 I V x - - 

2ojC xp 2 VJ C 2 xp , 



C xp = -jC xp + ojV p - ojV x - 2777 V x - 



Cxpt 

(A6) 



from which it can be verified that 



V x 

Vu 



c 



xp 



yOO 

v°° 



1/2, 
1/2, 
0, 



(A7) 



when t — y oo. From Eqs. (|A5I) and (IA7|) . we know that 
the state of the controlled harmonic oscillator tends to 



where \ij) e ) and \ip g ) are Gaussian states with first and 
second-order quadratures in the phase space as: 



j,e, (P)g,e — Pg.e 



) 2 



9,e' 



(P 2 ) 9 ,e - (P)le = V* e , 
(^) gte -{x) g , e {p) g , e =C*%. 

The first-order quadratures x g ^ e , p 9if , satisfy the equa- 
tions: 



dx„ 



d Pg,e = 



7 



X g . e dt + (A od =F X) Pg.edt 
1 

2 



^[Vle--] dW, 



dY, 



7 

(A od =f x) x gfi dt - -p g , e dt 

(ktY gte - k 3 Y g 3 e + k ) dt + y^C^dW, 

~ {Yg, e - X g . e ) dt - J~J=dW, 



1 



and the second-order quadratures V g e: Vg P e , C g f e satisfy 
the equations: 



yp 

g.e 



7 



- 1 Vg* e + 2(A od Tx)CZ + i 

1 

2 



-2^77 ( Vie ~ ~ n 



- 7 V^ e - 2 (A od T X) CZ + |-2t77 {C%Y , 
- 7 C^ e + (A od T X) V g p e - (A od T X) V g x e 
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Here, as an approximate estimation, we omit the higher- 
order disturbance induced by V ge , V?, and C g p , then 
the coefficients Pij(t) are given by (see, e.g., Ref. 36]): 

Pgg (t) = Pgg (0) , p ee (t) = Pee (0) , p ge (t) = p* eg (t) , 

where 72 is the dephasing rate of the qubit without mea- 
surement; and E(t), 0(i) can be expressed as: 

E W = X / [^e(s)p ff (s) -pe(s)Xg(s)]ds, (Bl) 



6(t) = X / [^(s)a; 3 ( S )+p e ( S )p 9 ( S )]d S . (B2) 
Jo 

By tracing out the degrees of freedom of the harmonic os- 
cillator, it can be shown that there exists a measurement- 



's, e 



induced dephasing factor exp [—£(£)]. 

As is proved in Appendix [A] we have 
Pg.e ~> Pg^e m the long-time limit. Thus, for any s > i > 
I/7, it can be approximately estimated as x gi e(s) ~ ^^e; 
Pg,e(s) ~ P^ e - Then, the measurement- induced dephas- 
ing after f can be approximately calculated as: 



CX P [%e(s)p g (s) - p e (s)x 9 (s)] d,S 

exp{[ x oor-^r)] (*-*)}> 



which leads to Eq. (p?Tj) . 



[1] H. M. Wiseman and G. J. Milburn, Quantum Measure- 
ment and Control (Cambridge University Press, Cam- 
bridge, 2009). 

[2] D. D'Alessandro, Introduction to Quantum Control and 
Dynamics (Chapman & Hall, Boca Raton, 2007). 

[3] H. Mabuchi and N. Khaneja, Int. J. Robust Nonlinear 
Control 15, 647 (2005). 

[4] D. Y. Dong and I. R. Petersen. larXiv:0910. 2350k 1 [quant- 
ph]. 

[5] N. Ganesan and T. J. Tarn, Phys. Rev. A 75, 032323 

(2007) ; R. B. Wu, T. J. Tarn, and C. W. Li, Phys. Rev. A 
73, 012719 (2006); D. Y. Dong, C. B. Zhang, H. Rabitz, 
A. Pechen, and T. J. Tarn, J. Chem. Phys. 129, 154103 

(2008) ; W. Cui, Z. R. Xi, and Y. Pan, Phys. Rev. A 77, 
032117 (2008); W. Cui, Z. R. Xi, and Y. Pan, J. Phys. A: 
Math. Theor. 42, 025303 (2009); M. Zhang, H. Y. Dai, 
Z. R. Xi, H. W. Xie, and D. W. Hu, Phys. Rev. A 76, 
042335 (2007). 

[6] J. Zhang, C. W. Li, R. B. Wu, T. J. Tarn, and X. S. Liu, 
J. Phys. A: Math. Gen. 38, 6587 (2005); J. Zhang, 
R. B. Wu, C. W. Li, T. J. Tarn, and J. W. Wu, Phys. 
Rev. A 75, 022324 (2007); J. Zhang, Y.-X. Liu, C. W. 
Li, T. J. Tarn, and F. Nori, Phys. Rev. A 79, 052308 

(2009) ; J. Zhang, R. B. Wu, C. W. Li, and T. J. Tarn, 
IEEE Trans. Automat. Contr. 55, 619 (2010). 

[7] V. P. Belavkin, J. Multivariate Anal. 42, 171 (1992); V. 

P. Belavkin, Commun. Math. Phys. 146, 611 (1992); V. 

P. Belavkin, Theor. Probab. Appl. 38, 573 (1993). 
[8] H. M. Wiseman and G. J. Milburn, Phys. Rev. Lett. 70, 

548 (1993); H. M. Wiseman and G. J. Milburn, Phys. 

Rev. A 47, 642 (1993). 
[9] A. C. Doherty, S. Habib, K. Jacobs, H. Mabuchi, and 

S. M. Tan, Phys. Rev. A 62, 012105 (2000). 
[10] M. R. James, H. I. Nurdin, and I. R. Petersen, IEEE 

Trans. Automat. Contr. 53, 1787 (2008); J. E. Cough, 

Phys. Rev. A 78, 052311 (2008); N. Yamamoto, Phys. 

Rev. A 74, 032107 (2006). 
[11] R. Ruskov and A. N. Korotkov, Phys. Rev. B 66, 

041401(R) (2002); L. Bouten, R. van Handel, and M. R. 

James, SIAM J. Control Optim. 46, 2199 (2007); L. 



Bouten, M. Guta, and H. Maassen, J. Phys. A: Math. 

Gen. 37, 3189 (2004). 
[12] H. Mabuchi and A. C. Doherty, Science 298, 1372 (2002). 
[13] R. L. Cook, P. J. Martin, and J. M. Geremia, Nature 

446, 774 (2007). 
[14] G. G. Gillett, R. B. Dalton, B. P. Lanyon, M. P. Almeida, 

M. Barbieri, G. J. Pryde, J. L. O'Brien, K. J. Resch, S. D. 

Bartlett, and A. G. White, Phys. Rev. Lett. 104, 080503 

(2010). 

[15] J. Q. You and F. Nori, Physics Today 58 (11), 42 (2005); 
Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. 
Phys. 73, 357 (2001); J. Clarke and F. K. Wilhelm, Na- 
ture 453, 1031 (2008); A. A. Clerk, M. H. Devoret, S. M. 
Girvin, F. Marquardt, and R. J. Schoelkopf, Rev. Mod. 
Phys. 82, 1155 (2010). 

[16] R. van Handel, J. K. Stockton, and H. Mabuchi, IEEE 
Trans. Automat. Contr. 50, 768 (2005); H. F. Hofmann, 
G. Mahler, and O. Hess, Phys. Rev. A 57, 4877 (1998); 
J. Wang and H. M. Wiseman, Phys. Rev. A 64, 063810 
(2001); X. S. Liu, W. Z. Liu, R. B. Wu, and G. L. Long, 
J. Opt. B: Quantum Semiclass. Opt. 7, 66 (2005); H. M. 
Wiseman, S. Mancini, and J. Wang, Phys. Rev. A, 66, 
013807 (2002). 

[17] S. Mancini and H. M. Wiseman, Phys. Rev. A, 75, 
012330 (2007); C. Hill and J. Ralph, Phys. Rev. A, 77, 
014305 (2008); A. R. R. Carvalho and J. J. Hope, Phys. 
Rev. A 76, 010301 (R) (2007); S. Mancini and J. Wang, 
Eur. Phys. J. D 32, 257 (2005). 

[18] L. K. Thomsen, S. Mancini, and H. M. Wiseman, Phys. 
Rev. A 65, 061801 (R) (2002); X. G. Wang, A. S. 
S0rensen, and K. M0lmer, Phys. Rev. A 64, 053815 
(2001); X. G. Wang, A. Miranowicz, Y.-X. Liu, C. P. 
Sun, and F. Nori, Phys. Rev. A 81, 022106 (2010). 

[19] P. Zhang, Y. D. Wang, and C. P. Sun, Phys. Rev. Lett. 
95, 097204 (2005); J. Q. You, Y. X. Liu, and F. Nori, 
Phys. Rev. Lett. 100, 047001 (2008); M. Grajcar, S. Ash- 
hab, J. R. Johansson, and F. Nori, Phys. Rev. B 78, 
035406 (2008); Y. D. Wang, Y. Li, F. Xue, C. Bruder, 
and K. Semba, Phys. Rev. B 80, 144508 (2009); F. Xue, 
Y. D. Wang, Y. X. Liu, and F. Nori, Phys. Rev. B 76, 



12 



205302 (2007); N. Lambert and F. Nori, Phys. Rev. B 

78, 214302 (2008); S. H. Ouyang, J. Q. You, and F. Nori, 

Phys. Rev. B 79, 075304 (2009). 
[20] R. Ruskov, K. Schwab, and A. N. Korotkov, Phys. Rev. 

B 71, 235407 (2005). 
[21] A. Hopkins, K. Jacobs, S. Habib, and K. C. Schwab, 

Phys. Rev. B 68, 235328 (2003). 
[22] J. Zhang, Y. X. Liu, and F. Nori, Phys. Rev. A 79, 

052102 (2009). 

[23] K. Jahne, C. Genes, K. Hammerer, M. Wallquist, E. S. 

Polzik, and P. Zoller, Phys. Rev. A 79, 063819 (2009); 

Ya. S. Greenberg, E. IPichev, and F. Nori, Phys. Rev. B 

80, 214423 (2009). 
[24] M. P. Blencowe and E. Buks, Phys. Rev. B 76, 014511 

(2007) ; E. Buks and M. P. Blencowe, Phys. Rev. B 74, 
174504 (2006); E. Buks, S. Zaitsev, E. Segev, B. Abdo, 
and M. P. Blencowe, Phys. Rev. E 76, 026217 (2007). 

[25] S. Etaki, M. Poot, I. Mahboob, K. Onomitsu, H. Ya- 
maguchi, and H. S. J. van der Zant, Nat. Phys. 4, 785 

(2008) . 

[26] T. Bhattacharya, S. Habib, and K. Jacobs, Phys. Rev. 

Lett. 85, 4852 (2000); S. Habib, K. Jacobs, and K. 

Shizume, Phys. Rev. Lett. 96, 010403 (2006). 
[27] M. A. Armen and H. Mabuchi, Phys. Rev. A 73, 063801 

(2006). 

[28] I. Siddiqi, R. Vijay, F. Pierre, C. M. Wilson, M. Met- 
calfe, C. Rigetti, L. Frunzio, and M. H. Devoret, Phys. 
Rev. Lett. 93, 207002 (2004); I. Siddiqi, R. Vijay, M. 
Metcalfe, E. Boaknin, L. Frunzio, R. J. Schoelkopf, and 
M. H. Devoret, Phys. Rev. B 73, 054510 (2006). 

[29] A. Lupascu, C. J. M. Verwijs, R. N. Schouten, C. J. P. M. 
Harmans, and J. E. Mooij, Phys. Rev. Lett. 93, 177006 
(2004); A. Lupascu, E. F. C. Driessen, L. Roschier, 
C. J. P. M. Harmans, and J. E. Mooij, Phys. Rev. Lett. 
96, 127003 (2006); A. Lupascu, S. Saito, T. Picot, P. C. 
de Groot, C. J. P. M. Harmans, and J. E. Mooij, Nat. 
Phys. 3, 119 (2007). 

[30] A. Isidori, Nonlinear Control Systems, 3rd edition 



(Springer Verlag, London, 1995). 
[31] K. Jacobs and A. P. Lund, Phys. Rev. Lett. 99, 020501 

(2007) . 

[32] S. Ashhab, J. Q. You, and F. Nori, Phys. Rev. A 79, 
032317 (2009); S. Ashhab, J. Q. You, and F. Nori, New 
J. Phys. 11, 083017 (2009); S. Ashhab, J. Q. You, F. 
Nori, Phys. Scr. T137, 014005 (2009). 

[33] J. K. Stockton, R. van Handel, and H. Mabuchi, Phys. 
Rev. A 70, 022106 (2004). 

[34] K. Petersen, Ergodic Theory ( Cambridge Studies in Ad- 
vanced Mathematics) (Cambridge University Press, Cam- 
bridge, 1990). 

[35] J. Combes, H. M. Wiseman, and K. Jacobs, Phys. Rev. 

Lett. 100, 160503 (2008). 
[36] J. Gambetta, A. Blais, M. Boissonneault, A. A. Houck, 

D. I. Schuster, and S. M. Girvin, Phys. Rev. A 77, 012112 

(2008) ; D. WahyuUtami and A. A. Clerk, Phys. Rev. A 
78, 042323 (2008). 

[37] A. C. Doherty and K. Jacobs, Phys. Rev. A 60, 2700 
(1999). 

[38] H. M. Wiseman and A. C. Doherty, Phys. Rev. Lett. 94, 
070405 (2005). 

[39] C. J. Hood, T. W. Lynn, A. C. Doherty, A. S. Parkins, 

and H. J. Kimble, Science 287, 1447 (2000). 
[40] C. J. Hood, M. S. Chapman, T. W. Lynn, and H. J. 

Kimble, Phys. Rev. Lett. 80, 4157 (1998). 
[41] T. Puppe, I. Schuster, A. Grothe, A. Kubanek, K. Murr, 

P. W. H. Pinkse, and G. Rempe, Phys. Rev. Lett. 99, 

013002 (2007). 

[42] A. Blais, R. S. Huang, A. Wallraff, S. M. Girvin, and R. J. 
Schoelkopf, Phys. Rev. A 69, 062320 (2004); A. Wallraff, 
D. 1. Schuster, A. Blais, L. Frunzio, R.-S. Huang, J. Ma- 
jer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, Nature 
431, 162 (2004). 

[43] J. Gambetta, A. Blais, D. I. Schuster, A. Wallraff, L. 
Frunzio, J. Mayer, M. H. Devoret, S. M. Girvin, and 
R. J. Schoelkopf, Phys. Rev. A 74, 042318 (2006). 



